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^ , Using Environmentally Friendly Renormalization, we present an analytic calculation of 

^ Y ' the series for the renormalization constants that describe the equation of state for the 

^ ^ ' 0{N) model in the whole critical region. The solution of the beta-function equation, for 

' the running coupling to order two loops, exhibits crossover between the strong coupling 

, fixed point, associated with the Goldstone modes, and the Wilson-Fisher fixed point. 

f**^ ■ The Wilson functions 7;^ , 7^ and 7^2 , and thus the effective critical exponents associated 

with renormalization of the transverse vertex functions, also exhibit non-trivial crossover 
between these fixed points. 
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1. Introduction 

The equation of state for the 0{N) model remains a subject of great interest (see, 
for instance, Refs. [TJ [2] for recent reviews). It exhibits crossover behavior between 
three distinct asymptotic regimes — the critical region approached along the criti- 
cal isotherm, the critical region approached along the critical isochore, and, finally, 
the coexistence curve. For iV = 1, the longitudinal correlation length remains finite 
away from the critical point on the coexistence cu rve, while for TV > 1, the existence 
of Goldstone bosons leads to infrared singularities .''^^ The problem of encapsulating 
these distinct scaling behaviors within one overall scaling function has been solved 
through an ab initio derivation from an underlying microscopic model.^ Specifically, 
in that work we obtained the equation of state for the 0{N) model using only the 
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Landau-Ginzburg-Wilson (LGW) Hamiltonian by implementing an Environmen- 
tally Friendly Renormalization (EFR) Group*^ which tracks the crossover between 
the fixed points that control the different asymptotic regimes. 

EFR is a formalism within a general framework of perturbative renormalization 
and the renormalization group, specifically designed to describe crossover phenom- 
ena, where the effective degrees of freedom at different scales can be quite distinct, 
leading to different scaling regimes and associated exponents. To do this, a suc- 
cessful renormalization should track the evolving nature of the effective degrees of 
freedom as a function of scale and that as the latter depend on the "environment" 
the reparametrization chosen should also depend of it. To illustrate this, if one con- 
sidered an interacting field theory in a three dimensional box of size L, one could 
renormalize the theory in an L independent fashion. When one considered physics 
on scales k, ^ one would find that the theory was perturbatively ill defined, 
whereas an appropriate L dependent renormalization made perturbative sense. The 
reason for this, of course, is that the effective degrees of freedom in the system are 
explicitly L dependent. An L independent renormalization ignores this important 
physical fact. The only fiuctuations being absorbed into the renormalized param- 
eters in this case are L independent, no matter what renormalization scale one 
chooses. L here is the parameter which induces the crossover and therefore a good 
renormalization scheme should be L dependent. In principle, basically any system 
will exhibit crossover behavior in some regime. Some pertinent examples are: sys- 
tems with a bicritical point^, bulk-surface crossoveri^ and dimensional crossover^. 

The advantage of this method relative to standard RG techniques is that it 
describes perturbatively the crossovers between any and all fixed points as opposed 
to the perturbative regime around one single fixed point. The disadvantage is that 
the Feynman diagrams that enter in the perturbative calculations are computed in 
the relevant environment where the finite part not just the asymptotic divergence 
is crucial. Applying EFR to study the crossovers inherent in the equation of state 
we obtain explicit functional forms that obey all required analyticity properties and 
require no phenomenological input, only the three Wilson functions 7a, 7<^, and 7^2 
which are deduced from within the theory. In particular, the case A'' = 1 was treated 
analytically in the one- loop approximation in Ref. [5] In a recent studjC^ we carried 
out the task again in the one- loop approximation, but now to include all A'^ > 1. 

In this paper we continue the extension towards the two-loop approximation 
of the universal equation of state by deriving the Wilson scaling functions to this 
order. The main motivation for doing this is to contribute with the crossover func- 
tions required to obtain the equation of state, and to calculate the first non-trivial 
correction to the critical effective exponent r], which follows directly from the Wil- 
son function Such a result is useful also as a checkout as no analytic form can 
be accessed to this order even for A^ = 1. As a considerable part of the formidable 
task of getting the equation of state for the 0{N) model to order two loops, we 
consider it worthwhile to present the calculation in this publication. To begin with. 
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m Section [2] we outline the renormalization group representation for the ab initio 
formulation that will be considered throughout this paper. Then, in Section |3] we 
obtain explicit expressions to the order of two loops for the bare vertex functions, 
the renormalization constants and the Wilson functions. In Section |4] we derive the 
two-loop order Wilson functions and solve the beta-function equation of the run- 
ning coupling, whose curve exhibits the crossover between the Wilson-Fisher and 
strong-coupling fixed points for A'^ > 1. The corresponding crossovers in the Wilson 
functions, and therefore in the effective critical exponents, are also presented in this 
section. Finally, we draw concluding remarks in Section [H 

2. Renormalization Group Representation 

The model is described by the standard LGW Hamiltonian with 0{N) symmetry 

n[v] = I d'^x Qv^'^V^'^-f ir(x)^V^+ , (1) 

which describes an A^-component scalar field ip m a d-dimensional space. Here r 
denotes the bare mass parameter, r = rc + ts, with rc being the value of r at 
the critical temperature Tc and Ib ~ ^21^21^^ where A is the microscopic scale. 
The value rc can naturally be interpreted in statistical mechanics as rc ccTc — Tm, 
where is the critical temperature predicted by mean field theory. As is well 
known, an additive renormalization for rc is first necessary to compensate for the 
critical temperature shift; then a further multiplicative renormalization oi ts is 
needed. 

The generator of connected correlation functions, W, is given by 

W[Ha]^lnZ, (2) 

where Z is the functional integral over the order parameter fields i^", with Hamil- 
tonian ([1]) and an external source Ha{x) 

Z[H]= [[dip]e-^^^^+^^"^"-^\ (3) 



The connected correlation functions are obtained by repeated functional differenti- 
ation of W with respect to ^^(x) and Ha{y). We denote these by 

Gil'^JJx^,...,XN,y,,...,yM). (4) 

In the same fashion, the vertex functions ri^'^ai are obtained by functional differ- 
entiation with respect to ip{x) and t{x) of the effective action r[i^], which is given 
as the Legendre transform 

rm = -W[Ha] + J d^X Ha{x)^''{x), (5) 

where (p"" (x) is the physical magnetization of the system defined as 
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In the ordered phase two types of modes exist: those along the external field Ha 
and those perpendicular to it. If we denote by n° to the unit vector in the direction 
of the external field, then by using the projectors 

p;'' = n'^n^ = S"'' - n°n^ (7) 

we can write a general vertex function as ^^g^fP i- When all subscripts are equal to 
£ or we compact them to one. For instance, v\'^'^ ^'^ is denoted rj^'^^\ Addition- 
ally, if there are no tp^ insertions (i.e. M = 0), the second superscript is omitted. 
That is, we write r(^'°) = F^^l 

As a consequence of the Ward identities of this model, all vertex functions can 
be expressed in terms of the transverse vertex functions. For instance, from the 
equation of state ri^^ = Ha we have Pj^^ — and V^p = H, so that use of the Ward 
identity F^^^ = F^^^ yields 

Ff)^ = iJ, fW=0. (8) 
Decomposing V)^^ produces F^ , FJ ' and F^ ^ . Ward identities then imply 

rf =FP) + Y^^ and Ff)=G. (9) 

Analogously, one may express any vertex function in terms of the y['^'^^\ In this 
sense, the transverse vertex functions are the building blocks of the theory. The 
Wilson functions 7i, to be defined shortly, can in particular be written in terms of 
these functions.!^ 



Due to the existence of large fluctuations in the critical regime, a renormalization 
of the microscopic bare parameters of the form 

t{mt,K) = Z^2{K)tBimt), (10) 

A(k) = Zx{k)Xb, (11) 

^{k) = Z-'/^k)^b, (12) 

must be imposed, where k is an arbitrary renormalization scale and rrit is the inverse 
transverse correlation length. The renormalized parameters satisfy the differential 
equations 

d 



— 

an 



= 7a(k)A(k), 



where 
where 
where 



7a(k) = 



K-r- In Zx 



J^{k) = K— \nZ^ 



(13) 
(14) 
(15) 



where on the right-hand side are the Wilson functions associated with this coor- 
dinate transformation and the derivatives are taken along an appropriately chosen 
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curve in the phase diagram, which we here denote by c. In this paper we are inter- 
ested precisely in finding the ji as crossover scaling functions to the order of two 
loops. 

Integration of the RG equation for any multiplicatively renormalizable r^^'^^ 
yields 

rr'^)(i, A,^) = e/r*(f 7.-M7..)^rr'^)(t(«), A(«),^(«)). (16) 

The renormalization constants Z^, Z^2 and Z\ are fixed by imposing the explic- 
itly magnetization-dependent normalization conditions on the transverse correlation 
functions 

dp2rf\p,t{K,K),X{K),ip{K),K) =1, (17) 

r^')(o,f(K,«),A(K),^(KO,«) = i, (18) 

r(^)(0,t(K,K),A(K),(^(«),K) =A, (19) 



while the condition 



=Tf\0,t{K,K),X{K),<p{K),K) (20) 



serves as a gauge fixing condition that relates the sliding renormalization scale k to 
the physical temperature t and magnetization (p. Physically, k is a fiducial value of 
the nonlinear scaling field m/ . 

Besides mj , the other nonlinear scaling field we use to parametrize our results is 



3 s„.r(^) 



(21) 



p2=o 



which is RG invariant. It represents the anisotropy in the masses of the longitudinal 
and transverse modes and is related to the stiffness constant ps = <fi'^ 9^2 

p2=0 

via = ^AyOg. With this renormalization prescription one may determine the 
Wilson scaling functions in terms of the nonlinear scaling fields rrit and the 
transverse and longitudinal propagators that appear in all perturbative diagrams 
can be parametrized in terms of them. 



3. Perturbative Series 

3.1. The bare correlation functions 

Within the ab initio formulation we are using, there appears the difficulty of calcu- 
lating Feynman diagrams that are no longer simple numbers but functions of the 
variables mj and instead. In the one-loop approximation it is possible to eval- 
uate the integrals involved analytically; however to higher orders this is no longer 
possible and this, numerically, complicates derivation of the equation of state. 
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In this section we present a perturbative expansion of the correlation functions 
and show expUcitly their dependence on the nonhnear scaling functions. We start 
from the two-loop effective action 



Tp\x) {A + r)Tp\x) + -T^{x) 



+ 41 ' "'^ 



2,G} {x, x) + 2{N -I) Ge {x, x) Gt {x, x) 
+ {N^ - 1) (x, x) ] - ^ y d^x d^y r {x) [SGI (x, y) 
+ {N-l)Ge{x,y)GUx,y)]v'{y), 



(22) 



which embodies the physics of the system. Notice that we are implicitly working 
with bare quantities and to simplify the writing we will be using diagrammatic 
notation. We represent the longitudinal propagator by a solid line and the transverse 
propagator by a dotted line: = p^+r+^lp'^ and G^^ — p^+r+^Tp^ respectively. 
It can be observed that in the symmetric ordered phase (magnetization ^ = 0) the 
two propagators are equivalent. Moreover, for models in the Ising universality class 
(TV = 1) the terms mixing propagators have no contribution in the correlation 
function. 

To obtain the two point correlation function, we take two functional derivatives 
of r[^] respect to the order parameter Tp. After a large but otherwise direct calcu- 
lation, we find the two-loop approximation for the two point correlation function in 
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momenta space 



r(2) („) = -± 



4! 



00 + 5<"-"0' 



-2(A'-l)^<Q, + -(«=-!) 



A3 
'4! 

4 



12 



^(iv-DOO+'^^-^^OO 



36 



- 18A 







e 



X2 



2*1 — ^Ajo — 2 A 



+0" 

+0")'''''-^(0-'€>°) 



I 



\ 277^2 



2 2 

3 ^ 3 



1 



+ ^0) I^V 



(23) 



where p denotes the external momenta flowing through the diagrams, the index 
indicating evaluation of the corresponding diagram at the point — 0. By using 

the projectors (l7|) onto this expression one may identify both the transverse Fi (p) 

(2) 

and the longitudinal F), {p) two point correlation functions. The transverse two 
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point vertex function is given by 



A^2 A 



6 ^ 



A 

'4! 
2 



' 36 



N 



1 



2 r 



- A(^2 



3 k::^ 



(24) 



where we have added the terms arising from the one-loop contribution. The deriva- 
tive of this function respect to the external momenta p is also required. For this we 
find 



9„.rp) (0) = 1 



36 



(iV-1) 



60 -9A^^0 



A2 



36 







(25) 



To obtain this result we have taken the derivative and then evaluated at the point 



P 



0. We emphasize this fact with a diagonal line crossing the diagrams. 



The four point correlation function can be calculated by taking four derivatives 
in Eq. ([22]) or alternatively by using the Ward identity By using the latter we 
get to 



.(4),nA_ ^ / _ N-1 



rr (o) = A-^A^( ^ 



A3 



12 



00 



+6 



-4 {N - 1) .Q. ^ + ^ (iv2 - 1) ,Q.2 



9 



3x3 

+ 2^ 



A3 
12 



4 /'"N 



3 v:^' 



-2A^ 



-2 ^2 



1 



2 



(26) 



(2,1) 



To complete this stage of the calculation we need the correlation function F 
which can be obtained if one derivative of the function F^ (p) respect to the param- 
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eter t is taken. The resulting expression is 



(2,i)_ A 



A2 
4! 



--(7V^-1) ( 2, 



" 4 



+ (iV - 1) 



A2 r 



18 



7^2 



— Xtp 



2^ 



(27) 



3.2. Wave function renormalization 

We now move on to consider the wave function renormahzation. To proceed, we 
shaU use the relationships that define both scaling variables and in the one- 
loop approximation, and then perform an inversion of these at the same order. The 
next step is to replace the resulting expressions into the one-loop terms within the 
correlation functions and then make the expansion up to two loops. By perform- 
ing these expansions into the first two terms of Eq. (j26p we get to the four point 
correlation function 



r(^) (0) = - |a| 



+ 3(^-l)< 



o 

2 
9 



N - 1 



X3 



+ A3 



\ —2 

-^B^B 



+ {N-1) 



X3 

12 



-2\Bip% 



+ f (^-1) 
A^- 1 ^ 



o 



I \ 4 TTT^ 

;> + 2 Ab'/'b 



O 



O 



(28) 
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(2 1) 

By repeating the same procedure into the one-loop terms of Fj ' , one obtains 



f ''^ (0) = 1 



A 



TV- 1 



2-10' 3 

2 



4! 



3(iV-l)<-><3> + -(JV>-l)< 

2, 



X2 



1 



+ {N-1) 



18 



2z: 



XBip% 



\ — 2 



+2^ 



+ ^(iV-l) 
A/'-l 



O 



+ 



3 

-o 

O 



(29) 



Using these expressions we shall calculate the renormalization constants that are 
necessary to find the Wilson functions to the order of two loops. 



3.3. Renormalization constants 

In this section we calculate the renormalization constants or Z functions. These are 
defined in terms of the bare correlation functions that we found in the previous 
section. For the function Z\ we can write 



Za = 1 - sAl 



o 



+ 



N -1 



X2 



+-(iV-l),Q,<^ + -(iV2-l) 



-Xb^b ( ^ + ^ 
9K 



+ (iV_i)^ 
^ M2 



3 '■■■h^^-'' 



-2Ab^|(^^ 



■^B 

' 18 



0-9Xb^10 



+ f (iv-i) 



2^ 



o 



O 



TV- 1 



O 



whereas for Z ^ the expression is 



36 



n0 



(7V-1) 



\2 

36 



(30) 



(31) 
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and similarly for Z^l one finds 



z-} 



Xb 
' 2 



o 



+ 



N -1 



4! 



+ |(7V-l),Q.^ + |(iV2-l) 



+ 



X2 



+ (iV-l) 



x2 r 
18 



1^ -2 



\ —2 

Xb<Pb 



+ 2, 



X2 



36 

-^BlfiB 



^0) 



2 r 



+ 



2^ 



O 



iV- 1 



o 



o 



(32) 



Notice that we have been working with the bare coupling As , which we now write 
explicitly in these expressions. 



3.4. The Wilson functions 

The Wilson functions are defined, in terms of the renormalization constants Ztp 
and Z^2 that we found in the previous section. From the definition of 7a, by using 
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Eq. (|5D|) . one obtains 



7A 



O 



N - 1 



+5<"-"00-i^<'^'-"< 



3^2 



\2 

+ (Ar_ i)_ 
^ '^12 



3 



3 

18 



1 



3^ 



\2 

^ '^18 
+ ^AVi?K 



o o 



2 

3^ 



A^2 



where we have used d/dhik := to simphfy the writing. 

Similarly, from the definition of and using Eq. pip , we get 



A^ r 
36 



6£) 



'0 



+ (JV-1) 



2 r 



36 



2D 



0, - A(^2^, 
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and finally from Eq. we obtain for the Wilson function 7^2 , 



7v= 



A ( ^ TV- 1 



1 



3 



i?. ( 2.^ + I - A^2^, 



A2 r 



36 



-9A^^^.0 



+ (TV - 1) 



A2 
36 



4(iV-l) 



iV- 1 _ 



o 



2 



4 » i ^ ^ 



(35) 



Notice that we have applied the I?k to the renormalization constants and then 
replaced A^ in terms of the renormalized coupling A. Also note that in these ex- 
pressions all the diagrams are given in terms of the nonlinear scaling fields nit and 
m^p explicitly. 



4. 3D Two-loop Results 

As anticipated in Section [U we are interested in obtaining expressions for the 7^ as 
crossover scaling functions. The natural variable is z, and so the next stage consists 
in writing the Wilson functions in terms of the nonlinear scaling fields z — mt/m^p. 
For instance, the Wilson function 7A which is cubic in the coupling (see Eq. (|33| ) 
turns out a quadratic in terms of z. It is no difficult to express the Feynman diagrams 
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appearing in the Wilson functions in terms of z. By defining the functions 

^■-1 ^ 
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3d \_ N-1 

87r2 I 9 



1 



-3/2^ 



/2(2) = -D, 



3 _ 2 1 



2O n^^-^)00-i8(^-^)< 
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(1 + ^) 
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1 f 2tt 



167r3 I 3 
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k{z) 
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6D 

3d 1 



2 



727r2 



1 + ^ 



11/, 1 
1 + ^ 



-3 



friz) = AV^« 



O 



+ 



A^- 1 



4 » 



2 



O 



3d 



2567r2 ^2 

(A^-1) 



-3/2 



1 + ^ 



-5/2 



/8(^) 



1 



A<^ 



the crossover scaHng functions 7a and 7,^ can be written respectively as 



7a {z) = -^/lA + (^/2 + ^/3 + ^^^^5^/4 + + ^^5^/6 

+ pr + {N-l)fs)x\ 



I'piz) 



(36) 

(37) 
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Further, definition of the functions 
/ TV - 1 



3d 1_ N-1 

"8772! 3 



3d 



r(i/2)^ 

(47r)3 
2(7V- 1) 



1 

" 18 
N ~ 1 



(TV -1)1?,/ 



9 



-1/2 



18 

iV- 1 



-3/2 



53(2) 



D 



3d 



(47r)3 



All I 1 
^ 1 + ^ 



12 / 1 
1 + ^ 



TT 1.4 

9+- 



54(2) = 



55(2) = 



~ o 



77^2 



allows one to write 



1^' 



in the form 



Iv^ (^) = - ^.9iA - 



1 iV- 1 1 
.92 + 7:33 H 7^54 + TT^/s 



18 



36'' 



^-1 , , Ay2 



(38) 



(39) 



Note that all the crossover scaling Wilson functions are finite in the limits of small 
and large values of z. For the functions (z) and gi [z) for which it is not possible 
to obtain an exact analytical expression, we have verified using Mathematica that 
their behavior is finite in both asymptotic limits. 

From the definition of 7a and its relation to the /3 function, we get the /3-function 
equation for the dimcnsionless coupling A: 



= -(4-d)A + 7AA. 



(40) 



As discussed in Ref. [SJ one attempts to reconstruct the Wilson functions from their 
series to a given order in the loop expansion using some method of resummation. 
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This is necessary as the direct perturbative series arising from the EFR scheme are 
divergent asymptotic series for any couphng strength, so that physical results can 
only be accessed from resummation procedures. Fade resummation is one technique, 
based on rational functions having the same power series expansion as the original 
series to the given order, that has been successfully used and is widely accepted for 
resumming perturbative series since the early works of Baker,^] and so we use it 
here to obtain a resummed series for the beta function; i.e. right-hand side of Eq. 
(|40| . A noteworthy fact is that using different Fade approximants one may estimate 
errors in the resummed series. This is especially useful for higher loop computations. 

There also exist more sophisticated techniques that can be used for resummation 
of asymptotic series, such as the Fade-Borel and conformal mapping methods. The 
former applies the Fade approximation to the Borel transform, whereas the latter 
is an improvement to the Fade-Borel method based on the mapping of a complex 
plane into the unit disk. More efficient techniques are based on re-expansions of the 
asymptotic truncated series in terms of special basis functions which are chosen to 
possess precisely the analytic behavior responsible for the divergence of the original 
series. These, and the method of variational perturbation theory, which is a system- 
atic extension of a variational approximation to path integrals, require information 
on the behavior of the series, so they are suitable for higher order expansions (see 
Ref. [12] for a review of these methods). 

For our two-loop calculations, we use for simplicity the Fade approximant. That 
is, we solve numerically the [2/1] Fade-resummed differential equation arising from 
Eq. (HUl) . The solutions we find are shown in Fig.[T] There, one may see the coupling 









250 
1 200 


X 








N=2 


















\ 150 










N=3 




\\ioo 










N=4 












In 














z 








-10 




-5 






5 


10 



Fig. 1. The Pade- resummed coupling A(z), for several values of A'^. Note the crossover to 
the strong fixed point for z — > 0. 

parameter interpolating continually between the fixed points of the model. The 
presence of the Wilson-Fisher fixed point A can be observed for large numerical 
values of z and, for > 1, the fixed point associated with the coexistence curve 
in the limit z -> 0. Notice the relative decrease of the maximum, on the left of the 
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curve, as N increases. More precise values in the asymptotic regimes, obtained from 
higher order approximations that introduce external numerical data, are known in 
the literature. This leads one to expect the same qualitative behavior of the curve 
to higher orders in the loop approximation. For completeness, in Table[l]we present 
the asymptotic values for the coupling in the limit z — >■ and in Table [2] the values 
corresponding to the Wilson-Fisher fixed point. However, as these values come from 



Table 1. Asymptotic values 
of the 3D coupling and Wilson 
functions in the limit z ^ 0. 
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Table 2. Asymptotic values of the 3D coupling 


Wilson func 


tions and critical exponents 


at the WF critical point. 




N A* 


7l 


7; 




/3 


5 


V 


7 


1 29.01 


1 


0.033 


0.43 


0.33 


4.8 


0.64 


1.26 


2 25.26 


1 


0.033 


0.52 


0.35 


4.8 


0.68 


1.33 


3 22.23 


1 


0.032 


0.59 


0.37 


4.8 


0.71 


1.39 


4 19.76 


1 


0.031 


0.64 


0.38 


4.8 


0.73 


1.45 



a two-loop calculation, no greater precision is expected. What is worthwhile noticing 
is the fact that, within our ah initio calculation and without external input data, 
we capture the crossover-function character of the running coupling interpolating 
from one fixed point to the other. By substituting the numerical solution for the 
coupling into the Pade-resummed 7a iz) function, we obtain the behavior of the 
Wilson functions showing the continuum crossover between the Wilson-Fisher fixed 
point and the fixed point associated with the coexistence curve. 

The crossover scaling function ^\(z) in Fig. [2] provides information on the effec- 
tive dimension of the system^. Given that the dimension employed to evaluate the 
expressions is c? = 3, asymptotically the value is 1, as expected. Nevertheless, it is 
of theoretical and experimental interest the local minimum values of this function 
for > 1. For iV = 1, the limit 7a shows mean field behavior as the fluc- 
tuations are suppressed. On the other hand, the crossover scaling function 7<p(z) 
directly corresponds to the effective exponent 77, see Fig. |3l As one would expect, 
in the Wilson-Fisher fixed point, this takes values that approximate those known 
from higher order calculations. We stress once again that, in the two-loop order 
approximation that we consider in this work, the contribution comes from the ah 
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Fig. 2. The Wilson function 7a (z) showing that 7^ — > 1 in the hmits z — > 00 and z — J- 0, 
except for A'^ = 1. 

initio calculation and no attempt has been made at improving numerical values us- 
ing external data. The numerical values that we obtain are, nevertheless, provided 
in table [5J Once more, in the limit 2; — s> the fluctuations are suppressed and this 
can be observed in the zero value of this anomalous exponent. Finally, the crossover 





1 ^ 


-10 ^"'^^ iV=T\-. 

\\\a'=4 , 

\ \f/N= 


/ -Q-Q2 
' -0.04 
-0.06 

3 

-0.08 



Fig. 3. The Wilson function 7(^(2). 



scaling function 7(^2 in Fig. |4] also shows non-trivial crossover. In the asymptotic 
limit z — 00, this function is related to the critical exponent v. It can be observed 
that the curves 71^2(2), for > 1, tend to follow the trajectory described by the 
Ising model in the limit z — s> 0, but in the end they separate from it. Again in this 
case, just as a reference, the asymptotic values of this function and for completeness 
the values of the critical exponents have also been included in both tables. 

5. Conclusions 

By using EFR it is possible, from an ab initio calculation, to derive a parametric 
form for the equation of state of the 0{N) model that has all desired analiticity 
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Fig. 4. The Wilson function 7^2(2). Notice the mean field behavior of the Ising model, 
AT = 1, in the limit z —> whereas 7^2 — >• 1 for iV > 1. 



properties in the three distinct asymptotic regimes. To order one loop it is even 
possible to obtain an analytic expression for iV = 1 in dimension three. We see that 
the fundamental building blocks for calculations are the Wilson functions. However, 
despite the fact that they play a privileged role they enter into expressions for physi- 
cal quantities in a non-trivial way, such as in integrals. In the case where the Wilson 
functions must be computed numerically in the first place this involves numerical 
subtleties that are not present in standard calculations, such as perturbation ex- 
pansions in e or l/N for problems where, for instance, the upper critical dimension 
or the symmetry of the order parameter change, respectively. In reality this is to be 
expected, calculating a crossover scaling function is much more complicated than 
calculating an exponent. 

From the EFR formalism it would seem that once the Wilson functions are cal- 
culated at a given order, the derivation of the equation of state at the same order 
should be somehow straightforward. However, the Feynman diagrams appearing 
in the Wilson functions to order two loops are crossover functions themselves and 
the difficulty precisely resides in systematically sum diagrams whose divergences do 
cancel. This is a problem we are currently sorting out. In this paper we have per- 
formed explicit two-loop order calculations of the transverse correlation functions, 
the renormalization constants and the Wilson functions. By solving numerically the 
beta-function equation, we have captured the crossover between the critical fixed 
point and the fixed point associated with the coexistence curve. The Wilson func- 
tions also show the crossover. In the limit z — ^ 00, the Wilson-Fisher fixed point is 
approached and 7i — >■ 7]^^ with 7a = 1 for d = 3. In contrast, in the limit z ^ the 
strong-coupling fixed point is approached and ji — >■ 7f*". For > 1 the Goldstone 
bosons dominate and 7a = 7,^2 = 1. For N = 1 however, this fixed point is mean 
field like as fluctuations are suppressed and 7, — > 0. 
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Appendix A. Derivatives of Feynman Integrals 

We summarize here expressions for some of the diagrams appearing in the numerical 
functions {z) and Qi (z). In terms of the functions 

f{x,y) = x{l- x) (1 -y) + y, 
g{x, y,z) = f {x, y) + y{l-x) z'"^, 

we can write 



Jo g{x,y,zy 



D^.^O:j:, = Aa I I dxdy- _^ 

Jo Jo g(x,y,z) 

where Aj, = - ^^^"^■'Jj^/"''^ K^ii-d> ■ Analogously, wc have 
Jo Jo g [x, y, z) 



{x{l-x)f""^{l~y)y^-''/^f{x,y) 



D^^—:^=Ba[ [ dxdy- _^ 
■v^-'- Jo Jo g{x,y,z) 

= ^ /' I' dxdy 
Jo Jo 



2 Jo Jo g{x,y,z) 



where Bj, = — ^^^j^ '^^ ^2(5-d) ■ Finally for the diagrams with the diagonal line 
crossing them, we find 



nl nl 

D^,^ = -Ad I J dxdy 
= -Bd ^ dxdy 
= -Bd I I dxdy 



y' 


-d/2 


{l-y){xil-x)f 


-'^'f{x,y) 






1 \5 — d 




x' 


-d/2 


{y{l-x)f-'^'{l- 


-y)f{x,y) 






g \x. y, z) 


1 


{x 


(1- 




-yf f{x,y) 



Jo g{x,y,z) 
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